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Nucleon sea structure functions are studied using Dokshitzer-Gribov-Lipatov-Altarelli- 
Parisi (DGLAP) equations with the massive gluon-quark splitting kernels for strange and 
charm quarks, the massless gluon-quark splitting kernels for up and down quarks, and 
the massless kernels for all other splitting parts. The SU (2) f fiavor symmetry for two 
light quarks, 'up' and 'down', is assumed. Gliick-Reya-Vogt(GRV) and Martin-Roberts- 
Stirling(MRS) sets are chosen to be the base structure functions at Qq = 3 GeV^. We 
evolve the sea structure functions from Qq = 3 GeV^ to = 50 GeV^ using the 
base structure function sets and DGLAP equations. Some (about 10%) enhancement is 
found in the strange quark distribution functions at low x(< 0.1) in leading order of the 
DGLAP equations compared to results direclty from those structure function sets at the 
the value of = 50 GeV'^. We provide the value of k and also show the behavior of 
k{x) = 2s{x) / {u(x) + d{x)) after the evolution of structure functions. 



The nucleon structure function at low x and at high x may bare non-perturbative 
effects. However, we can still investigate the structure function at not so low 
10~^ — 10"'^) with perturbative QCD. At this low x, sea quarks and glu- 
ons play important roles in QCD processes. Disagreements between theoretical and 
experimental results in various sum rules suggest the necessity of more improved 
theoretical analysis in nucleon structure functions, especially, in sea quark distribu- 
tion functions in this region of x. 

The traditional analysis of sea quark distribution functions for three light quarks 
(up, down, and strange) has been based either on the SU{3)f flavor symmetry or 

on an ad hoc choice of s{x) = ^u{x) = ^d{x) (or s{x) = u{x) = d{x)). Now, we are 
able to probe these sea structure functions at low x (up to ~ 10~^) and high 
with better statistics. The validity of the SU (3) / symmetry assumption and the ad 
hoc choice of the strange quark structure function has been at stake □. 

Although we believe that there are many signatures that treating strange quark 
on an equal footing as other two light quarks is not correct, it is not clear what 
physics should be applied and how to implement it to differentiate its roles in the 
structure functions. In this letter, we treat strange and charm quarks massive. The 
mass effect is one of several possible factors affecting the structure functions such 
as the non-perturbative QCD effect and the different Pauli exclusion effect for u, 
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d, and s quarks which should be taken into account due to different numbers of 
valence u and d quarks in the proton. On this-principle, however, there have been 
speculations that the effects may be marginal eI. 

To observe sea quarks in the nucleon (proton), we consider the process of the 
deep inelastic scattering of charged lepton off the nucleon. Howeve*, we do not con- 
sider the charged weak-current process for|-6implicity and clarity □. To obtain sea 
structure functions, we use massive kernels El for strange and charm quarks in gluon- 
quark^plitting kernels in Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) equa- 
tions □ while we keep the SU{2) f flavor symmetry for two light quarks, up and down, 
by taking massless kernels for the splitting function. We also keep massless kernels 
for all other splitting parts. Unlike our earlier work 13 where we used an approximate 
expression for DGLAP equations for the quark structure functions, we solve numer- 
ically and explicitly the DGLAP equations for the relevant structure functions in 
leading order of as- 

The DGLAP equations and the corresponding kernels to the quark and gluon 
distribution functions are given below. 
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Where, qi is the quark distribution function for the corresponding flavor 'i' and 'i' 
covers all possible quark and anti-quark flavors. We include only four quark flavors 
in this numerical calculation by excluding bottom and top quarks for simplicity. 
The inclusion of these quraks is straight forward but it will not change our current 
result noticably. 

We use massless kernels, Pg^, Pgg, Pgq in the calculation. For the gluon-quark 
splitting kernel, however, we use the massive kernel P™ for strange and charm 
quarks while we use the massless kernel Pqg for up and down quarks. Note that 
P™ Pqg as ^ (or m'^/Q'^ — > 0). Also, the terms in the massive kernel 
P™ have the power (m^/Q^) suppression like the ones appearing in higher twist 
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calculation. We choose the massive kernel, P™, for heavy flavors only in the g ^ q 

splitting part since g(x, Q^) is the most dominant structure functions at low x. The 
'+' prescription is defined in standard way, 



To solve numerically the DGLAP integro-differentioal equations for the quarks 
and gluons, we take the advantage of the differential form of the equations. We start 
the evolution of the structure functions at f = to and plug them into the original 
equations at next t — Iq + At. Then, iterate them until we get the desired value of 



As the Ip^Sf structure functions we take two typical structure function sets (GRV 
and MRS) EI, □ at the initial value, Qq and then, evolve the functions using the above 
DGLAP equations to the given value, Q^. The Gliick-Reya-Vogt(GRV) set contains 
structure functions in first two lowest orders in a^. On the other hand, the Martin- 
Roberts-Stirling(MRS) set contains the structure functions in higher orders in as 
but not in leading order. The structure functions in the leading order are used in 
this analysis. Therefore, for the same order of the calculation, it is more consistent 
to use the GRV set for the evolution and compare the evolved result with the 
corresponding GRV structure functions. However, we still include the result using 
the MRS set for comparision. 

We start the evolution at the value of Qq = 3 GeV^ and end it at = 50 
GeV^. There is no special reason to choose the value of Qo = ^ GeV^. This choice 
of the Qq value is only to get the non-zero charm distribution function from the 
GRV set in the perturlsative region. 

Our goal is to find the low x behavior of heavy sea quark (strange and charm) 
structure funcitons. We set strange and charm quark masses to be rus — 0.2, 
0.5 GeV and nic — 1.45 GeV. We also include the evolution result (at Q^ = 50 
GeV^) with the strange quark mass iris = GeV for the comparision to the strange 
structure functions of GRV and MRS at the same Q^. We choose the value of Aj^ 
to be Aj^ — 0.2 GeV for Nf — 4 which agrees with the value of Aj^ used in the 
GRV structure function set while it is not much different from the one in the MRS 
set in which Aj^ = 0.23 GeV. 

The expression of v in the massive kernel in eq. (||) contains the mass threshold 
(for pair creation process) at parton level. In other words, > implies that 
Q2(1 -z)/z> 4m^ like that the Bjorken variable x = Q^/[Q^ + (W^ - M^)] should 
be restricted by the mass threshold, = Q'^{l~x)/x > '^m^ for the heavy quark 
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ruQ where M is the negligible mass of the nucleon and W is the invariant mass of 
the virtual photon and the nucleon. 

We observe about 10% enhancement in the strange quark structure function at 
low X and the given . The boost of this distribution function at such low x is due 
to the introduction of the massive kernel in the DGLAP equations. This massive 
kernel effect persists for the surprisingly high value of for the strange quark mass 
iris. 

Due to the complicated dependence of rrig/Q^ and z, P™ is not so suppressed as 
was originally thought (See figures 1 and 2). However, the effect from the massive 
kernel can not be strong for the charm quark (See figure 3). It is fine to use the 
massive kernel for all range of that we considered, — 3 GeV^ ^ 50 GeV^. 

To show this mass effect more closely, we also use the massless kernel for the 
strange quark. As shown in fig. 1, the result of the evolution with this massless 
kernel is smallest and thus closer to the strange quark structure function of the 
lowest order GRV structure function set at = 50 GeV^. 

Similar behavior can be found in using the MRS structure functions set (in fig. 
2). The evolved strange structure function with the zero strange mass is smallest in 
low X compared to the strange structure functions with the non-zero strange mass. 
The difference is that all the obtained strange structure functions after the evolution 
are smaller than the one directly from the MRS set. Remind that the MRS set we 
used is in next to leading order in as while our result is in leading order in as like 
the GRV set used here. Therefore, the next to leading order correction may be 
bigger than the mass effect at low x. 

We applied the same method to the charm quark structure function (in fig. 3). 
We obtained the evolution result which is smaller in using both the GRV set and the 
MRS set. The reason for this behavior is that the mass threshold of the charm quark 
pair(cc), 2.9 GeV, is big enough to suppress the corresponding structure functions. 

The ratio, k{x) = 2s{x) / {u{x) + d{x)) , is also important in judging the degree of 
the SU (3) / symmetry breaking in terms of variable x. Note the different definition 
of k(x) from the usual one k, 

^ 2(5(^)) ^ '^lLo^^(^)dx 
i{u{x)) + {d{x))) J^^^ 

xu{x)dx + J^^Qxd{x)dx 

Here, the integration range should be [0,1] instead of [x,l]. However, due to the 
unknown nature of the structure functions at low x, we have to use a low value of 
X close to '0' but not x = 0. 

k{x) contains some uncertainty at low x like the structure functions do. However, 
this k{x) contains uncertainty also in high value of x, x > 0.5, in our calculation, 
due to very small and thus uncertain values of the sea structure functions in the 
range of x. For example, all sea structure functions, s{x), u{x), and d{x), are too 
small with big uncertainty for x > 0.5 for k to be meaningful. As shown in the 
figure 4, k(x) at the high value of x is not reliable. Fortunately, in calculating k, 
however, those structure functions at this range of x give the negligible effect. This 
uncertainty comes from the technicality in calculation and is nothing to do with the 
nature of QCD unlike the previous uncertainty. 

We found interesting behavior in k(x). k{x) stays almost constant in x after 
the use of the MRS set while it changes somewhat (it shoots up fast at low x) after 
the use of the GRV set as shown in the figure 4. The value of k also shows the 
interesting feature. Values of k both from GRV and MRS sets stay around 0.4 as 
shown in the table (|l|) and increase slowly for the smaller integration range of x 
(i.e., bigger value of Xmax in the table). The value of Xmax should be '1'. Remember 
the definition of k in eq. (||). We chosed the lower limit of x in the integration to 
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be '0.001'. On the other hand, k's from our calculation stay around '0.6'. This is 
a big difference. The experimentally measured number of n may bring the way to 
improve our analysis further. 

Table 1. First column, Xmax, stands for the maximum value of the integration. For example, 
0.1 value of x corresponds to the integration range of x of (0.001, 0.1) and 0.5 to the range 
of (0.001, 0.5). Second, third, fourth, and fifth columns in the table stand for the values of k 
using the structure functions directly from the GRV set, evolved results with the GRV functions, 
the structure functions directly from the MRS set, and evolved results with the MRS functions, 
respectively. In all cases, = 50GeV^ . 
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Figure 1: The strange structure functions arc shown here with different strange 
masses after the evolution with the GRV base structure function set. The upper 
sohd, large dashed, small dashed, dotted, and lower solid curves correspond to the 
evolved strange structure functions with strange mass, rus = .5GcV, TOs = .2GeV, 
nis = OGeV, and the GRV strange structure functions with = 50GeV^ and 
= 3GeV^ respectively. 
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Figure 2: The strange structure functions arc shown here with different strange 
masses after the evolution with the MRS base structure function set. The upper 
sohd, long dashed, short dashed, dotted, and lower solid curves correspond to the 
evolved strange structure functions with strange mass, nis = .5GeV, TOs = .2GeV, 
nis = OGeV, and the MRS strange structure functions with = 50GeV^ and 
= 3GeV^ respectively. 
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Figure 3: The charm structure functions are shown here with the mass nic = 
1.45GeV for the GRV and the MRS base structure function sets. The upper sohd, 
upper long dashed, short dashed, dotted, lower solid, and lower long dashed curves 
correspond to the GRV charm structure functions with — 50GcV^, the MRS 
charm structure functions with = 50GeV^, the evolved charm structure with 
Q2 _ 50GeV^ and the GRV set, the evolved charm structure with = 50GeV^ 
and the MRS set, the GRV charm structure functions with = 3GeV^, and the 
MRS charm structure functions with = SGeV^ respectively. 
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K(x) 



Kappa from eGRV at Q'^2=50 GeV'^2 
Kappa from eMRS at Q'^2=50 GeV^Z 
Kappa from GRV at 0^^2=50 GeV'^2 
Kappa from MRS at Q'^2=50 GeV'^2 
Kappa from GRV at Q'^2=3 GeV'^2' 
Kappa from MRS at 0^^2=3 GeV'^2 
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Figure 4: k{x) in terms of x with two structure function sets and the evolved struc- 
ture functions with two different sets. The soHd, long dashed, short dashed, dotted, 
long dot-dashed, and short dot-dashed curves correspond to n{x) from the evolved 
structure functions with the GRV set and = 50GeV^, the evolved structure fmic- 
tions with the MRS set and = 50GeV^, the relevant structure functions of the 
GRV set with Q"^ = 50GeV^, the relevant structure functions of the MRS set with 
(32 = SOGeV^, the relevant structure fimctions of the GRV set with = 3GeV^, 
and the relevant structure functions of the MRS set with = SGeV^. 



